o 
o 



- 1—1 

X 



Minimum Enclosing Area Triangle 
with a Fixed Angle 

Prosenjit Bose* Jean-Lou De Carufel^ 
September 16, 2010 



Abstract 



o 

<D 

Given a set S of n points in the plane and a fixed angle < uj < it, we show how to 
find all triangles of minimum area with angle uj that enclose S in 0{n log n) time. We 
also demonstrate that in general, the solution cannot be written without cubic root. 

o 

U 1 Introduction 

i In geometric optimization, the goal is often to find an optimal object or optimal placement of 

an object subject to a number of geometric constraints (see |AS98j for a survey). Examples 
include finding the smallest circle enclosing a point set [Meg83, Wel91] or finding the smallest 
circle enclosing at least k points of a point set of n points {k < n) |ESZ94l IMat95j . In our 
setting, we study the following problem: given a set S of n points in the plane, find all 
the triangles of minimum area with a fixed angle uj, < u < tt, that enclose S. When 
no constraint is put on the angles, Klee and Laskowski |KL85j gave an 0(n log 2 n) time 
algorithm for finding the minimum area enclosing triangle. This was later improved to 
0(n log n) by O'Rourke et al. [OAMB86J which is optimal. Bose et al. [BMSSedJ provided 
j> optimal algorithms for the setting where one wishes to find the minimum area isosceles 

triangles. The setting we explore here is in between the two. We place a restriction on the 
angle but do not insist on the triangle to be isosceles. Our solution, which we outline below, 
uses ideas from the solutions of Klee and Laskowski and Bose et al. 

The five main steps of the algorithm are presented in Sections [2j [3j |3J [5] and [6j Each 
section outlines one step, proves the mathematical formulas involved and gives its time 
complexity. As we present in Section [5j at some point, the algorithm needs to calculate 
the roots of a fourth degree polynomial. This is unfortunate when it comes to numerical 
robustness. However, in Section [7| we show that we cannot dodge such algebraic expressions. 
We finally conclude and present some related open problems in Section |8l 
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^Ottawa-Carleton Institute for Computer Science. 



1 



2 Overview and Preliminaries 



Since the solution to the general problem only needs to consider the vertices of the convex 
hull of S, the first step is to compute the convex hull. In the remainder of the paper, we 
assume that the input is a convex n-gon with vertices given in clockwise order. 

Let P be a convex n-gon (refer to Figure [I]) . We denote the edges and the vertices of P 




Figure 1: In this example, P is a quadrilateral and w = |tt. 

in clockwise order by e$ and V{ for < i < n — 1 (all index manipulation is modulo n). Here 
and in the following sections, as we present the algorithm, we trace each step through the 
example of Figure [TJ 

We begin with two definitions. 

Definition 1 (w-wedge). Let q be a point in the plane and u> be an angle (0 < ui < n). 
Let A and A' be two rays emanating from q such that the angle between A and A' is u. 
We say that the closed set formed by q, A, A' and the points between A and A' creates an 
u-wedge, denoted W(u, q, A, A'). The point q is called the apex of the uj-wedge. An uj-wedge 
W encloses a polygon P when P (1W and both A and A' are tangent to P. 

Definition 2 (w-cloud). Let P be a convex n-gon and u be an angle (0 < u < ir). By 
rotating an enclosing u-wedge around P , the apex traces a sequence of circular arcs that we 
call an w-cloud (refer to Figure^. 

There are many technical details involved in the solution to this problem. Before getting 
too caught up in the technical details, let us first review the general approach to our solution. 

Since we only consider enclosing triangles with an angle of u, each optimal triangle can be 
constructed from an w-wedge that encloses P. Therefore, we consider all possible cu-wedges 
that enclose P. The apices of these w-wedges lie on an w-cloud Q which consists of a linear 
number of pieces of circular arcs (refer to [BMSSedJ). Then, for each of these u- wedges, we 
find the minimal triangle by identifying a third side. For the triangle to be optimal, the 
midpoint of this third side has to touch P (see Proposition [3] and Corollary [6]) . Hence, for 
each enclosing o;-wedge, there is one and only one triangle to consider for optimality. 

Moreover, when the apex q of an cu-wedge W moves clockwise around the w-cloud fl, 
the midpoint m of the third side of the optimal triangle moves clockwise around P (see 
Lemma [7]). As W moves, we note the positions of q where m leaves an edge of P and where 
m enters a new edge of P. These positions of q are important event points. They divide Q 
into a linear number of sections (see Section El). Let S be one of these sections. We prove 
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that the minimum triangle having a vertex (hence an angle u) on S can be computed in 
constant time (see Lemmas [To| [TTj [l2]and 13). We then have a linear number of candidates 
(one for each section of fl) to consider rather than infinitely many if we had to take all 
possible enclosing a;- wedges into account. What remains is to identify the optimal ones from 
these linear candidates. With this in mind, the technical details will fall into place. 



Step 1 Compute the w-cloud around P and denote it by Q. 

Q consists of n! = 0(n) circular arcs that we denote in clockwise order by Tj for < i < n' — l. 
The intersection point of Tj and T i+1 is denoted by u i+ \ for < i < n' — 1 (refer to Figure [2]). 
We also refer to Tj by the closed set [ui, Wj+i] containing all the points of IY Step 1 takes 




Figure 2: Step 1: Q is the i7r-cloud of P {n 1 = 6). 



0{n) (refer to jBMSSedj ). 



3 Optimal Solution Given a Fixed uj- Wedge 

In this section, we present a routine that computes the minimum enclosing area triangle 
from a fixed u- wedge enclosing P (see Algorithm [4]). Then we describe Step 2 that uses 
this routine. 

Take any w-wedge W = W(u, q, A, A') that encloses P. Find b G A (respectively c 6 A') 
such that P is enclosed in Aqbc and the midpoint m of be is on P. We claim that 

1. it is always possible to find b and c satisfying these properties, 

2. Aqbc is the minimum triangle enclosing P that can be constructed with W(u, q, A, A'), 

3. it takes 0{n) time to compute Aqbc. 

These claims are proven in Proposition |3j Algorithm |4| Lemma [5] and Corollary [6] 

Proposition [3] addresses the case where q is outside of P. It was first proven by Klee and 
Laskowski (refer to [KL85J). We propose a slightly different proof. Among other things, we 
stress the fact that Aqbc always exists. 
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Proposition 3. Let P be a convex n-gon and q be a point outside of P. Let Aqbc be 
the triangle that circumb scribes P such that the midpoint m of segment be lies on P. The 
following is true: 

1. Aqbc exists and is unique (hence it is well-defined). 

2. Among all triangles that circumbscribe P and have q as a vertex, Aqbc is the unique 
one of minimum area. 

Proof: 

1. Take b and c on the rays that support P from q such that an edge e of P lies on be (refer 
to Figure |3]). If m G e, then we are done. If not, suppose without loss of generality 



q 




Figure 3: Aqbc exists. 



that m is between e and c. Move b along its ray such that b gets farther from q and 
be stays tangent to P. Therefore c gets closer to q. Since m moves continuously as 
b moves along its ray, m will eventually touch one of the P's edges. This continuity 
argument implies the existence of Aqbc. 

Suppose there are two triangles, namely Aqbc and Aqb'c' with midpoints m and m' 
respectively. We show that these two triangles are equal. There are three cases to 
consider: (1) m = m', (2) m 7^ m! and both are on the same edge of P or (3) m 7^ vn! 
and both are on different edges of P. 

(1) If m = ml is not a vertex of one of the edges of P, then it lies on the interior of 
one of the edges of P (refer to Figure [1]). Hence this aforementioned edge is in 
be H b'd . This implies that b — b' and c = d by construction. 

Suppose m = mf is the vertex of one of the edges of P. Triangles Ambb' and 
Amcd are congruent by the following: (refer to Figure [5]). 

• \b'm\ = \dm\ since m is the midpoint. 

• Zb'mb = Zdmc since they are vertical angles. 

• I bm I = I cm I since m is the midpoint. 
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edge of P 

Figure 4: m — m! is not a vertex of P. 



g 




Figure 5: m — ml is a vertex of P. 



Therefore Zmbb' = Zmcd and Zmb'b = Zmc'c. Hence we have the following: 

a = Zb'mb = Zc'mc , 
(3 = Zmbb' = Zmcc' , 
7 = Zmb'b = Zmc'c , 

where a+/3+7 = n. So Zbcq = a+7 because Zbcq and Zc'cb are supplementaries. 
Therefore Zbqc = TT — a — f3— , -f = 0, which is impossible unless a = 0. We conclude 
that b = b' and c = d . 

(2) Suppose m m' and they are on the same edge of P. If neither m nor w! is a 
vertex, then this situation is similar to the one in Figure |4| hence b — b' ', c — d 
and m = ml . This is a contradiction so this situation is impossible. 
Without loss of generality, suppose that ml is a vertex of P and that \m'b\ < \m'c\. 
Let b" be the point on the line through b'd such that Segments bb" and cd are 
parallel, (refer to Figure [6]). Triangles Am'bb" and Am'cd are similar by the 
following: 

• Zbm'b" = Zcm'd since they are vertical angles. 

• Zm'bb" = Zmlcd since they are alternate angles. 

Therefore, since \mlb\ < \m'c\, then \m'b"\ < \m'd\. However, since ml is 6'c"s 
midpoint, \m!b"\ > \m'd\. This is a contradiction so this situation is impossible. 
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edge of P 
Figure 6: m! is a vertex of P. 



(3) Suppose m m' are not on the same edge of P. This case is similar to the case 
where m 7^ ml and both points are on the same edge of P. 

2. Let Aqb'c' be a triangle of minimum area that circumbscribes P and has q as a vertex 
(refer to Figure [7]). Suppose for a contradiction that the midpoint m! of b'd does not 

3 




Figure 7: Aqbc has minimum area. 



lie on one of the edges of P (refer to Figure [7]) . 

Since P is convex and Aqb'c' circumbscribes P, the following construction is possible. 
Let v G b'd be the vertex of P closest to m! . Without loss of generality, suppose 
vd < vb'. Let b"d' be a line segment that goes through v and such that b" G qb' and 
d' lies on the line through q and c'. Let 9 = Zdvd'. Zb'vb" = Zdvd' since they are 
vertical angles. If 9 is sufficently small, then vd' < vb" and Aqb"d' circumbscribes P. 

Therefore 

area(Aqb'd) = area(Ab'vb") + area([qb"vd}) 

= -\vb'\\vb"\ sin(#) + area([qb"vc]) 

> -\vc'\\vc"\ sin(0) + area([qb"vc]) 

= area(Ac'vd') + area([qb"vd]) 
= area(Aqb"c") 
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which contradicts the fact that Aqb'c 1 has minimum area. So m' lies on one of the edges 
of P and Aqb'c' is a local minimum among triangles circumbscribing P and having q 
as a vertex. We proved in Point [T] that there exists one and only one such triangle. □ 

Given the setting of Proposition [3j we show how to compute Aqbc in 0(n) time. In 
Algorithm [4j the edges of P are considered in clockwise order. The variable edge indicates 
the edge of P that is currently being considered. Algorithm [4] uses Lemma [5] that is presented 
after this routine. 

Algorithm 4. (Minimum Enclosing Triangle with a Fixed u-Wedge) 

• INPUT: An n-gon P and an enclosing u-wedge W = W(w, q, A, A') that supports P 
at vertices vq and Vk (refer to FigureUft). 




Figure 8: Computing Aqbc takes 0(n) time. 



• OUTPUT: The minimum triangle enclosing P that can be constructed with W(u>, q, A, 
A'). 

1. edge < — e&. 

2. If the line through edge does not intersect A or A' , 

• edge < — next[edge). 

• Go to H 

3. Let b (respectively c) be the intersection point of the line through edge and A (respec- 
tively A'). If the midpoint of be is on P, return Aqbc. 

4- If the midpoint of be is to the left of edge, 

• edge < — next{edge). 

• Go to 
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5. If the computation reaches this step, it means that the midpoint of be is to the right of 
edge. 

• Place the cartesian coordinate system on W such that q = (0,0) and A is the 
positive x-axis (refer to Figure^. 

• Let v = (s, t) := edge PI previous(edge) . 

• Takeb = (2(s — tcot(ui)),0) and c = (2tcot(u),2t) (refer to the proof of Lemma^ . 

• Return Aqbc. 

In the worst case, Algorithm [i] considers all the edges of P. Since it spends 0(1) time 
per edge, it takes 0{n) time total. 

Lemma 5. Let A and A' be two lines intersecting at q and let u (0 < u < n) be the angle 
between A and A'. Let v be a point within W(o>, q, A, A') (v A and v £ A'). In 0(1) 
time, we can compute b G A and c G A' such that b, v and c lie on a single line and v is the 
midpoint of be. 

Proof: Without loss of generality, suppose A is the x-axis, A' is the line^a; — cot(u)y = 
and let v = (s,t) (t > and s — cot(u)t ^ 0) (refer to Figure [9]). 




Figure 9: Proof of Lemma 5 



Let A" be a line containing v such that A is not parallel to A" and A' is not parallel to 
A". The line A" is the one containing the points b and c we are looking for. There are two 
cases to consider: (1) A" is not vertical or (2) A" is vertical. 

(1) Suppose A" is not vertical. Let b be the intersection point of A and A", and c the 
intersection point of A' and A". The general equation of A" is y = X(x — s) + t 

1 In what follows, several algebraic expressions are written using cot(-) where it seems that they could be 
simplified by using tan(-) instead. However, we write cot(-) in order to properly deal with the angle |-7r. 
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for A 7^ and y 7^ cot (a;). Calculating the intersection point of A and A", and of 
A' and A", we find that the general coordinates of b and c are b = (^y^,0) and 

c = (^ 1 _^ C Qt( a; ) co ^( a; ); i-Acot(o,) )" ^ e are l°°ki n g for A such that \bv\ = \vc\, so we 
need to isolate A in 

Xs-t A 2 , n , 2 / t-Xs . , V f t-Xs 

which solves to A = 2fCQ ^- ) _ s - Therefore 6 = (2(s — £cot(a;)), 0) and c = (2£cot(a;), 2£). 
Note that A 7^ because t 7^ and | 7^ cot (a;) because s — cot(o;)£ 7^ 0. 

(2) Suppose A" is vertical. Therefore < uj < |7r, otherwise A" is not the expected line. 
Using the notation of the previous case, we get that the equation of A" is x = s. 
Moreover, b = (s, 0) and c = (s, stan(o;)). We are looking for \bv\ = \vc\, which means 

(s - s) 2 + (0 - tf = (s - s) 2 + (stan(w) - t) 2 . 

So this situation occurs if and only if 2t = stan(w) and the solution is b = (s, 0) = 
(2(s — tcot(u)), 0) and c = (s, stan(a;)) = (2tcot(cu), 2t). 

Hence the global solution is b = (2(s — tcot(co)), 0) and c = (2£cot(o;), 2t) in all cases. □ 
Note that what we obtained is more general than our claims. Because q is outside of P, 
Proposition [3] together with Algorithm [4] and Lemma [5] not only compute the minimum trian- 
gle enclosing P that can be constructed with W(uj, q, A, A'), it also computes the minimum 
triangle enclosing P and having q as a vertex. 

The case where q is on P is similar. We conclude with it. 

Corollary 6. Consider an u-wedge W = W(to, q, A, A') enclosing P such that q is on P. 
The minimum triangle enclosing P that can be constructed with W can be computed in 0(n) 
time. 

Therefore, given any fixed w-wedge W enclosing P, we can compute the minimum en- 
closing area triangle that can be constructed with W. The midpoint m of the third side 
of this optimal triangle has to be on P. In Step 3 (see Section we need to know the 
exact position of m for one fixed w-wedge enclosing P. However, it does not matter for what 
enclosing w-wedge we know the position of m, as soon as we know for one. Therefore, in 
Step 2, we fix an arbitrary u- wedge W enclosing P and we compute the position of m. 

Step 2 Let q G r be such that q = u . Consider the w-wedge W = W(cu, q, A, A') that 



encloses P (refer to Figure 10). Apply Algorithm E] with P and W. 



Consider Figure [2j Any w-wedge enclosing P and having its vertex in To \ {uq,u\} is 
such that Vq € A' and Vi 6 A. This property easily generalizes to the case where the apex 
of the w-wedge is on uq or on u\. Therefore, when we specify q & T and q = uq, we imply 
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6/ A 




Figure 10: Step 2: Aqbc is constructed with W(|7T, q, A, A') where g = uq € IV 



that t>o ^ A' an d v\ G A, and hence the w-wedge is unique. However, in the example of 
Figure 10, there are infinitely many other u- wedges enclosing P and having q = uq as a 
vertex to consider because uq is on a vertex of P (indeed uq = vq). The algorithm will 
consider these other to- wedges later (see Section EJ. 

From Proposition |3j Algorithm |4j Lemma [5] and Corollary [6j Step 2 takes 0(n) time. 



4 Walking Around the cj-Cloud 

In the previous section, we computed the minimum area triangle enclosing P that can be 
constructed with a given w-wedge W(o>, q, A, A'). The next step is to consider all possible 
w-wedges enclosing P and having their apex on Q. For each of these w-wedges, we compute 
the minimum triangle enclosing P. The final solution is the minimum among all these 
triangles. Of course, there are infinitely many such triangles. However, we show that we 
need to consider only 0(n) of these triangles. 

In what follows, we consider the points of Q and the related o;-wedges in clockwise order. 
When thinking of all possible w-wedges enclosing P and having their apex on Q, one must 
pay attention to the following case. If q is the endpoint of one of the circular arcs Ui of Q, 
and if at the same time this endpoint is one of the vertices Vj of P, then there are infinitely 



many w-wedges to consider (refer to Figure 11 ). This situation occurs when oj is greater than 



the angle between the two edges and ej. Then the u- wedge turns around the corner 
created by ej-i and ej while its apex q stays on the vertex Vj. 

The following lemma describes the behaviour of m (the midpoint of be constructed in 
Step 2) as q moves clockwise around Q. Some hypotheses of this lemma are written using 
the expression "close enough". Its meaning is twofold. First, it is easy to see that as an 
enclosing w-wedge W = W(oj, q, A, A') moves continuously such that its apex q stays on Q, 
the midpoint m of be also moves continuously. Secondly, let e be an edge of P. If m G e, it 
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Figure 11: |7r-wedge turning clockwise around a corner. 

is possible to move q such that m stays on e. Such a displacement of q has to be of less than 
e for an e > and such an e always exists. 

What needs to be shown is that m moves clockwise around Q. 

Lemma 7. As an u-wedge W = W(u, q, A, A') moves clockwise such that its apex q stays 
on Q, the midpoint m of be — the third edge of the minimum triangle enclosing P that can be 
constructed with W — moves clockwise along P. Specifically, take q', q" G Q with q" clockwise 
from q' . 

1. If q' 6 Tj \ {ui + \} (0 < i < n' — 1), q" <E]q',u i+ i], and q' and q" are close enough to 
each other, then we have the following: 

(a) If m! G Cj \ {vj,Vj + i} (0 < j < n — 1), then m" G]m',t> J+1 ] C ej. 

(b) If m! = v j (0 < j < n — 1), then m" G Cj (possibly m" = m'). 

2. If q' = q" = Ui (0 < i < n' - 1) and u t = v t (0 < I < n - 1 ). Let W and W" be two 
different co-wedges enclosing P with q' = q" = Ui as an apex. If W" is clockwise from 



W (refer to Figure 11) and W and W" are close enough to each other, then we have 
the following: 

(a) If m' G Cj \ {vj,Vj + i} (0 < j < n — 1), then m" & ]m',Vj + i\ C ej. 

(b) If m' = Vj (0 < j < n — 1), then m" G Cj (possibly m" = m'). 



PROOF: Without loss of generality, b'c' is on the (refer to Figure [L2]) . 



1. (a) If q' and q" are close enough to each other, then m" G ej. Therefore b" = (b",0), 
b' = (b' x ,0), c" = «,0) and d = (c' x .,0) with b" x < b' x < c" x < c' x . Hence m" = 

b''+c'' . b'+c' i a _ i / i ^_ 

x 2 x < x 2 x = m , so m G \m , fj+ij C ej. 

(b) If q' and q" are close enough to each other, then the only situation to discard is 
m" G ej_i \ {vj}. So suppose m" G e^\ \ {vj} for a contradiction. Therefore m' 
and m" both belong to ej_i Hence an argument similar to the one of the previous 
case leads to m" = Vj, which is a contradiction. 
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Figure 12: Proof of Lcmma[7] 



2. The proof is similar to the one of Point [TJ 



□ 



Therefore, the midpoint m of the third side of the triangle is either a vertex of P or on 
an edge of P. The goal of Step 3 is to identify the sections of Q where the midpoint is a 
vertex and the sections of Q where the midpoint is on an edge of P. 

Step 3 Move the w-wedge with apex q clockwise along Q and maintain b, c and m as defined 
in Step 2 (refer to Section [3]). Collect all of the following three types of event points 



(refer to Figure 13). 



q 3 = u 2 :m 3 = v 3 
q 2 = ui : m 2 = v 3 



<? 4 : e 3 = [m 4 ,«o] C 6 4 c 4 
q 5 =u 3 : m 5 = v 

= it4 : "16 = 



Qi ■ e 2 = [v2,mi] C b\C\ 
q = u :m £ e 2 
qi2 ■ e 2 = [mi2,v 3 ] C b 12 Ci 2 
qu : ei = [wi,mu] C 6nCii 
q w : ei = [mi ,w 2 ] Q bi c w 




v 



q 7 : e = [m7,ui] C 6 7 c 7 



q s =u 5 : m 8 e e 



g 9 : e = [^0,^9] C 6 9 c 9 



Figure 13: Step 3: qa = uq = vq is such that mo G e 2 . Note that uq = vo gives birth to two different event 
points, namely q^ and q\ 2 . They correspond to the ^7r- wedge enclosing P and having uq £ T as a vertex, 
and the |7r- wedge enclosing P and having wo £ T$ as a vertex. q 2 — u\ : m 2 = v 3 means that q 2 = iti, 
7742 = ^3 an d b 2 c 2 n P — {v 3 }. It is an event point of the first type. At such a place, m does not move even 
though q does. (74 : e 3 — [7714, Vo] C 64 C4 means that 7714 = ?J3 and e3 C 64 C4. It is an event point of the 
second type, qg : eo = [t>o,77ig] C bgCg means that my = v± and eo C &9C9. It is an event point of the third 
type. 
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Type 1 q is on the intersection point of two consecutive circular arcs of Vt for an i 
with < i < n' — 1. Formally, q = Ui for an i with < i < n' — 1. 

Type 2 g is such that the third side 6c of the triangle is on an edge of P (for an 
i with < i < n — 1) and the midpoint m of 6c is on the first vertex Vi of 
(when the vertices of P are considered in clockwise order). Formally, q is such 
that m = Vi and = [i>i, fj+i] = [m, fj+i] C 6c for an i with < % < n — 1. 

Type 3 q is such that the third side be of the triangle is on an edge e; of P (for an 
i with < i < n — 1) and the midpoint m of 6c is on the last vertex v j+i of 
(when the vertices of P are considered in clockwise order). Formally, q is such 
that m = v i+ i and = [i>j, = [v^ m] C 6c for an z with < i < n — 1. 

For each event point q iy save its type together with the location of mj, the midpoint of 
the third side of the triangle. 

It is easy to find the event points of the first type and the following two lemmas show 
how to find the event points of the second and third type. Lemma [8] helps identifying event 
points related to Case [T] of Lemma [7| As for Lemma [9j it helps identifying event points 
related to Case |2] of Lemma [7l 

Lemma 8. Let T = [v,w] be an arc of a circle, A" be a line and Vj G A" be a point. It is 
possible to find a triangle Aqbc such that q e T, b,c G A" and Vj is the midpoint of be (or 



decide that there is no such point) in 0(1) time (refer to Figure 




Figure 14: Illustration of Lemma |8] 



PROOF: Without loss of generality, T is the locus of the point q such that Zvqw = u. Hence 
we can take v = (0,0), and w = (2r sin(a;), 0), where r is the radius of T. Let a = ZvjVW, 



(3 = Zvjwv and q be any point of T (refer to Figure 15). There are three cases to consider: 



(1) either A" has strictly negative slope, (2) A" has non-negative slope (3) or A" is vertical. 

(1) Suppose A" has strictly negative slope. The general equation for A" is y = Xx + /i 
(A < 0) and hence Vj = (s, As + /i) for sGi Note 6 = ZqvjW. Also let A (respectively 



13 



Figure 15: Proof of Lemma 8 



A') be the line through q and v (respectively through q and w). Let b (respectively c) 
be the intersection point of A and A" (respectively of A' and A"). 

Trivially < 9 < n — uj, but there are other restrictions, c must be strictly to the right 
of Vj, so (3 — uj < 6. Also, b must be strictly to the left of Vj, so 9 < n — a. Finally, 
q must not reach the point where the line through q and w is parallel to A". Hence 
9 < arctan(A) — uj. To summarize, we have 



By elementary trigonometry, we get 

q = (2r cos(#) sin(6 l + u), 2r sin(6 l ) sin(6 l + u>)) . 

Then, finding the general equation of the line through q and v, and of the line through 
q and w, we can calculate the general coordinates of b and c: 



1 - Acot(0) ' 1 - Acot(0) J ' 
2r sin 2 (o;)cot(6 l ) + 2rsin(o;) cos(u;) + /i cos(u;)cot(6 l ) — fism(uj) 
sin (u) cot (9) + cos(u;) — A cos(w)cot(#) + Asin(o;) 
(2rsin(o;)A + /i)(sin(co')cot(6 l ) + cos(w)) \ 
sin(u;)cot(#) + cos(o;) — A cos(ci;)cot(^) + Asin(a;) J 



Therefore we want to find 9 such that v j is the midpoint of be, which leads to the 
equation 



< 9 < 7T - uj , 

f3 — uj < 9 < min(7r — a, arctan(A) — uj) . 



(1) 
(2) 



b 




1 / //cot(0) 

2 Vl - Xcot(9) 



+ 



2r sin 2 (o;)cot(6') + 2r sin(a;) cos(w) + /i cos(co')cot(6 l ) — (ism(uj) 



) 



sin(o;)cot(^) + cos(a;) — A cos(o;)cot(6') + Asin(co>) 



= s 
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or 

(2r sin 2 (c<;)A + 2/icos(w)A — (ism(uj) — 2sAsin(w) + 2sA 2 cos(o;))cot 2 (^) 
+ (2r sin(cu) cos(w)A — 2rsin 2 (o;) — 2ficos(u) — 2/isin(c<;)A + 2ssin(c<j) 

— 4s cos(w)A — 2sA 2 sin(a;))cot(^) 
+ sin(w) — 2r sin(u;) cos(w) + 2s cos(w) + 2sA sin(a;)) = . (3) 

Since this is a quadratic equation in cot(#), it is solvable in constant time for 9 satisfying 
([!]) and ([2]) (or it is possible to decide that there is no such solution in constant time). 

(2) Suppose A" has non-negative slope. Using the notation of the previous case and a 
similar argument, A > and 

< 9 < 7T - uj , (4) 
max(/3 — uj, arctan(A)) < 9 < n — a . (5) 

We obtain the same quadratic equation ^ that is solvable in constant time for 9 
satisfying Q and (JsJ) (or it is possible to decide that there is no such solution in 
constant time). 

(3) Suppose A" is vertical. Therefore the general equation for A" is x = s for s G KL 
< uj < \tt otherwise there is no solution. There are three subcases to consider: (3.1) 
s < 0, (3.2) s G]0,2rsin(u;)[ (3.3) or s > 2rsin(w). 

(3.1) Suppose A" is vertical and s < 0. Let Vj = (s,t) for a t < 0. Trivially < 9 < 
71 — u>, but there are other restrictions. Whit the notation of the previous cases, 
w must be between q and c, therefore 9 < \tx — uj. b must be strictly above Vj so 
9 < tt — a. c must be strictly under Vj, so (3 — uj < 9. It all sums up to 

< 9 , (6) 
(3 — uj < 9 < min — a, -n — uj^j . (7) 

By elementary trigonometry, we get 

q = (2r 008(6*) sin^ + uj), 2rsin(9) sin(9 + uj)) . 

Then, finding the general equation of the line through q and v, and of the line 
through q and w, we can calculate the general coordinates of b and c: 

s 

s, 



cot6» 

(s — 2r sin(w))(sin(a;)cot(^) + cos(cu)) 
cot(6') cos(w) — sin(w) 
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Therefore we want to find 9 such that Vj is the midpoint of be, which leads to the 
equation 

If s (s — 2r sin(u;))(sin(u;)cot(#) + cosfw))' ^ 



or 



2 \cot(0) cot(0) cos(w) — sin (a;) 



(2rsin 2 (u;) + tcos(w) — s sin(w))cot 2 (#) 
+(2r sin (a;) cos(a;) — 2s cos(cj) — t sin(c<;))cot(6>) 
+ssin(w) = . (8) 

Since this is a quadratic equation in cot(6 l ), it is solvable in constant time for 9 
satisfying ^ and (|7| (or it is possible to decide that there is no such solution in 
constant time). 

(3.2) Suppose A" is vertical and s G ]0, 2r sin(cu)[. In this case there is no solution. 

(3.3) Suppose A" is vertical and s > 2rsin(o;). Using the notation of the case "A" is 
vertical and s < 0" and a similar argument, 

max ( (3 — oj, -n J < 9 < n — a , (9) 



2 



< ir-u . (10) 



We obtain the same quadratic equation (|8j) that is solvable in constant time for 
9 satisfying (|9| and ( 10 ) (or it is possible to decide that there is no such solution 
in constant time). □ 

Lemma 9. Let q be a point, A" be a line and Vj G A" be a point. It is possible to find the 
triangle Aqbc such that Zbqc = uj, b,c G A" and Vj is the midpoint of be in 0(1) time (refer 
to Figure 16). 




Figure 16: Illustration of Lemma |9] 
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PROOF: Without loss of generality, A" is the x-axis, Vj = (0,0), q = (s,t) with t > and b 
is to the left of c. Denote by d the center of the circumcircle of Aqbc and by r, its radius. 

Suppose oj 7^ 2 71 "- By elementary geometry, d lies on the line segment bisector of be and 
Zbdvj = Zcdvj = uj (refer to Figure 17). Hence d = (0,h) for an h G K. Therefore, the 



?=(s,t) 




Figure 17: Proof of Lemma [9] 



equation of the circumcircle i^j 



Since g is on this circle, then 



x 2 + (y - /i) 5 



s 2 + (t - /i) 



/i 2 



cos^w) 



so 



COS w 



sin^w) 



cos^w) 

Ucos(w) - ^ s 2 sin 2 (w) + t 2 ^j 



from which 



teotfu;) — \l s 2 sin 2 (w) + t 2 csc(w), 



tcot(w) + y s 2 sin (lo) + t 2 csc(u;), 



If u = =7T, a similar argument leads to the same solution. 



□ 



There are n' = 0(n) event points of the first type. There are at most n event points of 
the second type because there are n edges. For the same reason, there are at most n event 
points of the third type. Therefore, Step 3 (computing the event points) takes 0(n) time. 



2 We need the assumption lu 7^ \-k here, otherwise co ^ u ^ is undefined 
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5 Finding the Optimal Solution when the Apex is on 
a Circular Arc 



Between any two consecutive event points, there is a single arc of a circle (that might be 
reduced to a single point if qi is on one of the vertices of P). Consider such a circular arc. 
As q moves along it, the minimum triangle enclosing P changes. In this section, we show 
how to compute the optimal triangle for a fixed arc of a circle between two event points. 
Four different cases can occur: 

• either the circular arc is reduced to a single point or not. If the former is true, then 
the apex q is forced to stay on one of the vertices of P. 

• Either the midpoint m is forced to stay on a vertex of P or not. 

The following lemmas describe how to compute the minimum triangle enclosing P in these 
four different situations. 



Lemma 10 describes how to compute the minimum enclosing triangle when q moves along 
an arc of a circle and m is forced to stay on one of the vertices of P. 

Lemma 10. Let F = [v, w] be an arc of a circle and Vj be a point. It is possible to find the 
triangle Aqbc of minimum area such that q ET, q, v and b lie on the same line, q, w and c 



lie on the same line, and Vj G be in 0(1) time (refer to Figure 18). 

q 




Figure 18: Illustration of Lemma 10 



PROOF: The strategy is to first fix a point q on T and then follow the proof of Lemma [5] 
in order to construct the minimum triangle for this fixed q. Then we move q along T while 
maintaining the minimum triangle. The smallest one among all these minimal triangles is 
optimal. 

Without loss of generality, T is the locus of the point q such that Zvqw = u. Hence we 
can take v = (0,0) and w = (2r sin(o;), 0) where r is the radius of T. Let q be any point 
on T (refer to Figure 19). Note 7 = \vjv\, 5 = \vjiv\, a = ZvjViv (0 < a < ir), (3 = Zvjiuv 
(0 < (3 < n and < a + (3 < ir) and 9 = Zqvw. Let w' be the orthogonal projection of w 
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Figure 19: Proof of Lemma 
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onto the line through q and v. Let v'j be the orthogonal projection of Vj onto the line through 
w and w'. Let A" be the tangent to T at q. Finally, note s = dist(vj, A") and t = \vjVj\. 
By elementary trigonometry, we get 

q = (2r cos(6 l ) sin(6 l + u), 2r sin(6 l ) sin(# + oj)) . 

By elementary trigonometry and geometry, we get 

s = \vjVj\ + \w'q\ 

= 5 sin ^-7r — 9 + (3^j + 2r sin(ci;) sin(0)cot(a;) 

= S cos(# - (3) + 2r cos(w) sin(6>) , 
t — 7 sin(7r — — a) 
= 7sin(^ + a;) . 

With respect to the proof of Lemma |5j let b (respectively c) be the point on the line qv 
(respectively on the line qw) such that \qb\ = 2(s — teot(o>)) (respectively |gc| = 2tcsc(a;)). 
Therefore Aqbc is of minimum area for this fixed q (and therefore for this fixed 9) by Lemma[5] 

We need to study how the area of Aqbc evolves as q moves along T. Trivially < 9 < 
7i — u, but there are other restrictions. Since Vj must not be on the line through qv, 9 < n — a. 
Vj must not be on the line through qw neither, so tt — (3 < 9. Finally, we must have \qv\ < \qb\ 
and \qw\ < \qc\, so arccot Q(cot(/3) — cot(a))) —u<9 and 9 < arccot (|(cot(/3) — cot(a))). 
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It all sums up to 



max ( 0,arccot ( -(cot(/3) — cot (a)) 



u < e 



9 < min yn — oj, arccot ^-(cot(/3) — cot (a)) 

(3 — oj < 9 < n — a . 

By the proof of Lemma [5j the area of Aqbc as a function of 9 is 

a(9) = 2(s - tcot(w))t 

= 2((5cos(9 -0) + 2r cos(w) sin(0)) - (7sin(# + a)) cot (a;)) (7 sin(# + a)) 
= 27 sin(0 + a) (5 cos(6 l — /3) + 2r cos(a;) sin(0) — 7cot(w) sin(# + a)) . 

With the sine law applied on Aw jW, one gets 

2r sin (a;) sin(/3) 



(11) 

(12) 
(13) 



7 
5 



sin(a + (3) 
2r sin (a;) sin(a) 
sin (a + (3) 



Therefore, 



a(9) = 27sin(6> + a)(5cos(6> - /3) + 2r cos(w) sin(6>) - 7cot(w) sin(6> + a)) 
8r 2 sin(a) sin(/3) sin (a;) . 



sin 2 (a + (3) 



sin(/3 — a; — 9) sin(a + 



We have 



a'{9) 



8r 2 sin (a) sin(/3) sin(c<;) 



sin (a + f3) 
8r 2 sin(a) sin(/3) sin(a;) . 



(sin(/3 — a; — 0) cos(a + 9) — cos(/3 — oj — 9) sin(a + 9)) 



sin(a - /3 + oj + 26>) . 



sin^a + /3) 

In order to find the candidates for minimum and maximum of a. we need to solve 



sin(a - (3 + oj + 29) = 



(14) 



for 9 satisfying ([11], (12) and (|13j). 

We look for solutions to a — (3 + oj + 29 = kn for k G Z. At first sight, we have 
^ _ kn-a+p-uj ^ we neec j ^ study this general solution in order for it to satisfy (11), (12) 
and (13). Putting this solution against these constraints, we get 

kir-a + /3-uj 

(15) 
(16) 
(17) 



max ( 0, arccot ( -(cot(/3) — cot (a)) ) — oj ) < 



kn — a + (3 — u 



< min ( 7r — oj, arccot ( -(cot(/3) — cot(a)) 



kn — a + (3 — oj 
[3 — oj < < 7r — a 
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which simplify to 



max ( ± ( — kir + a — f3 + 2 arccot I — (cot(/3) — cot(a)) 

u; < min (kir — a + f3, —(k — 2)tt + a — (3) , 
max (— kir + a + (3,(k — 2)tt + a + /3) < w . 



< w 



(18) 

(19) 
(20) 



If k < 0, (20) and (19) lead to — kn + a + (3<tu<kn — a + (3 < —kn + a + (3, which is a 



contradiction. If k > 2, (20) and (19) lead to (k - 2)tt + a + (3 < u < —(k - 2)vr + a - (3 < 



(k — 2)tc + a + /3, which is a contradiction. If k — 1, then 9 
provided that 



7T — a+f}— UJ 



is a valid solution 



max ( ± ( 7T — (i -i ./ — 2arcco1 ( -(cot(/3) — cot(a)) 



< a; < min(7r — a + (3,tt + a — (3) 



In this case, 9 



■K—a+/3—ui 



is a maximum as we show below. 



and 



a"{9) 



16r 2 sin(a) sin(/3) sin(o;) 



sin (a + (3) 



cos(a - /3 + u + 26>) 



ir — a + f3 — cu 



16r 2 sin(a) sin(,9) sin(a;) cos(7r) 
sin 2 (a + (3) 
16r 2 sin(a) sin(/3) sin(a;) 



sin (a + (3) 



< 



7T — Q + — W 



is a maximum. 



so 

Otherwise, (14) has no solution, therefore a(9) is monotonic. The conclusion is that the 
minimum of a is at one (or both) of the extremities of the domain of a. □ 



Looking at the proof of Lemma 10, it suggests that the minimum triangle might not 



exist. Indeed, if < (3 — u and n — a < tt — u, then f3 — oj < 9 < n — a by (11), (12) and (13) 



But since the minimum of a is at one (or both) of the extremities of the domain of a, then 
it does not exist. However, in the setting of the general problem, this does not occur. By 
Proposition |3j the minimum area triangle always exists for a given u- wedge so 9 will always 
vary inside an interval that includes its extremities. 



Lemma 11 describes how to compute the minimum enclosing triangle when q moves along 
an arc of a circle and m is forced to stay on one of the edges of P. 

Lemma 11. Let V = [v,w] be an arc of a circle and A" be a line. It is possible to find the 
point g 6 T such that the line through qv, the line through qw and A" form a triangle of 



minimum area in 0(1) time (refer to Figure 20). 
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A" 

Figure 20: Illustration of Lemma 11 




PROOF: Without loss of generality, T is the locus of the point q such that Zvqw = u. Hence 
we can take v = (0, 0), and w = (2r sin(cu), 0) where r is the radius of T. Let q be any point 



on T (refer to Figure 21). There are three cases to consider: (1) A" has strictly negative 



slope, (2) A" has non-negative slope or (3) A" is vertical. 

(1) Suppose A" has strictly negative slope. The general equation for A" is y = \x + // 
(A < 0). Note 9 = Zqvw. Also note A (respectively A') the line through q and v 
(respectively through q and w). Note b (respectively c) the intersection point of A and 
A" (respectively of A' and A"). 

Trivially < 9 < it — u, but there is one more restriction, v must be between q and b, 
hence 9 < arctan(A) — u. It all sums up to 

< 9 , (21) 
9 < n -co , (22) 
9 < arctan(A) - u . (23) 

By elementary trigonometry, we get 

q = (2r cos(6 l ) sin(6 l + u), 2r sin(6 l ) sin(6 l + u j) 

cos(o;)cot(6') + sin(cu)cot 2 (^) cos(u;) + sin(a;)cot(#) 

ZT* ZT 

cot 2 (#) + 1 ' cot 2 (#) + 1 

Then, finding the general equation of the line through q and v, and of the line through 
q and w, we can calculate the general coordinates of b and c: 

/icot(9) /i \ 

1 - Acot(0)' 1 - Acot(0) J ' 

2r sin 2 (a;)cot(0) + 2rsin(o;) cos(cj) + /i cos(w)cot(6') — /isin(o;) 



sin (a; ) cot (6* ) + cos(a;) — A cos(cu)cot(0) + Asin(w) 
(2rsin(w)A + yu)(sin(u;)cot(0) + cos(cu)) 
sin(a;)cot(0) + cos(a;) — A cos(cu)cot(0) + Asin(a;) 

Let X{6) = cot(0), therefore 

/ sin(a;)X 2 + cos(u;)X sin(a;)X + cos(a; 
Q = 2r v9 , ; ,2r 



b 



X 2 + l ' X 2 + 1 

—fxX —fi 
XX — 1' XX — 1 

(2r sin 2 (w) + /i cos(cj))X + 2r sin(co>) cos(c<;) — /i sin(aj) 
(sin(w) — Acos(a;))X + cos(w) + Asin(a;) 
(2rsin(a;)A + /j)(sm(u)X + cos(cj)) 
(sin(w) — A cos(a;))X + cos(w) + A sin(co>) 
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Here is the area a of Aqbc. 

a(X) 
= ^\Qb\\(lc\sm(uj) 

1 , . / ((2Ar sm(u) + /i)X 2 + 2r(A cos(w) - sm(u))X + a - 2r cos(w)) 2 
= 2 Sm(w) \/ (AX-1) 2 (1+X 2 ) X 

' ((2Ar sin(w) + /i)X 2 + 2r (A cos(w) - sin(w))X + /i - 2r cos(w)) 2 
((A cos(w) - sin(u;))X - A sin(w) - cos(w)) 2 (l + X 2 ) 

1 . ((2Arsin(o;) + /i)X 2 + 2r(Acos(w) - sm(u))X + fi - 2rcos(o;)) 2 

2 Sm ^' ) | (AX - l)((Acos(w) - an(u))X - Asin(w) - cos(w))| (1 + X 2 ) 

1 . ((2Arsin(o;) + /z)X 2 + 2r(Acos(w) - sin(w))X + // - 2rcos(w)) 2 

2 Sm ^' ) (1 - AX)(Asin(u;) + cos(u;) - (Acos(u;) - sin(w))X))(l + X 2 ) 



The reason why the absolute value | • | disappeared is twofold. Firstly, by (23), 9 < 
arctan(A) — uj < arctan(A) = arccot (y), so AX = Acot(6>) < 1. Secondly, by (23), 
9 < arctan(A) — u>, therefore 



9 < arctan(A) 
9 + uj < arctan(A) 



uj 



9 + uj < arccot 



A 



COt(6 +L0) 

cot (uj) cot (9) - 1 
cot(0) + cot(w) 
Acot(w)cot(#) — A 



> 



> 



1 

A 
1 

A 



< 9 + to < vr by (21) and (22) 



< 1 



cot(6>) + cot (a;' 

Acot(u;)cot(0) - A < cot(6) + cot (a;) 

Acot(w)X - A < X + cot (a;) 

(Acot(w) - 1)X < A + cot (a;) 

(Acos(cu) — sm(u))X < Asin(a;) + cos(a;) 



A < by the hypothesis, 



< 9 + to < 7r by (21) and (22), 



so cot(0) + cot(cj) > 0, 



< uj < n by the hypothesis, 
so sin (a;) > 0. 



Now we need to find for what values of X is cr(X) minimum, wich means that we need 
to find for what X does cr'(X) = 0. 



a'{X) 
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= — - sin(w) ((2Ar sin(w) + /i)X 2 + 2r(A cos(w) — sin(w))X + /i — 2r cos(cj)) x 

((—ft sin (a;) + 2r sin 2 (a;)A + 2A 3 r sin 2 (a;) + 2Xfi cos(u;) + fi\ 2 sin (a;) 
+ 4A 3 r cos 2 (u) - 4rA 2 cos(w) sin(u;))X 4 

+ (— 2rsin 2 (u;) — 2/icos(u;) — 6A 3 r sin(a;) cos(u;) + 2A 2 rsin 2 (a;) 

+ 2ft\ 2 cos(w) - 4/iA sin(w) - 12A 2 r cos 2 (w) + lOAr sinfw) cos(w))X 3 
+ 6r(cos(u;) + A sin(a;))(— sin(u;) + 2A cos(a;) + A 2 sin(a;))X 2 
+ (2rsin 2 (a;) — 4r cos 2 (a;) — 2yUcos(a;) — 10A 2 r sin 2 (u;) — 4yU,Asin(a;) 

+ 2ftX 2 cos(u;) + 2A 3 r sin(a;) cos(u;) — 14Ar sin (a;) cos(cu))X 
+ (At sin 2 (a;) A — ftX 2 sin(w) — 2Xfi cos(a;) + 2r sin(a;) cos(a;) 

— 2rA 2 cos(a;) sin (a;) + fism(u))) -r- 
((1 - AX) 2 (-A sin(w) - cos(a;) + XX cos(a;) - X sin(w)) 2 (l + X 2 ) 2 ) 

Therefore ft'(X) = if and only if 

(2Ar sin(a;) + ft)X 2 + 2r(A cos(w) - ain(u))X + fi - 2r cos(w) = (24) 

or 

(—ft sin (a;) + 2r sin 2 (a;) A + 2A 3 r sin 2 (a;) + 2Xfi cos(u;) + /aA 2 sin(a;) 

+ 4A 3 r cos 2 (w) - 4rA 2 cos(w) sin(w))X 4 
+ (— 2r sin 2 (a;) — 2ficos(co) — 6A 3 rsin(a;) cos(a;) + 2A 2 rsin 2 (a;) 

+ 2ftX 2 cos(u) - 4/xA sin(a;) - 12A 2 r cos 2 (a;) + lOAr sin(w) cos(w))X 3 
+6r(cos(co>) + A sin(w))(— sin(a;) + 2A cos(a;) + A 2 sin(a;))X 2 
+ (2r sin 2 (a;) — 4rcos 2 (co>) — 2/zcos(co>) — 10A 2 r sin 2 (w) — 4/iAsin(a;) 

+ 2/iA 2 cos(a;) + 2A 3 r sin (a;) cos(u;) — 14Ar sin(a;) cos(u))X 
+ (4r sin 2 (a;)A — ftX 2 sin(w) — 2Xfi cos(a;) + 2r sin(a;) cos(a;) 

— 2rA 2 cos(a;) sin(w) + fism(u))) = . (25) 

The first equation is quadratic in X = cot(#), therefore it can be solved in constant 
time for X and then it remains to solve for 9. The second equation is quartic in 
X = cot(#), therefore it can be solvec^jin constant time for X and then it remains to 
solve for 9. 

Actually, what we are really looking for is fi'(X(9)) = 0. Since 

fi'(X(9)) = fi'(X)X\9) = fi'(X)csc(9)cot(9) , 

9 = ^tt is also a candidate. And of course, the extremities of the domain of 9 are also 
candidates. 

3 By Galois theory, polynomials of degree d < 4 can be solved exactly in constant time. Refer to |DF03| 
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Overall, we get at most nine candidates for minimum (at most two from the quadratic 
equation, at most four from the quartic equation, 9 = and the extremities.) So it 
can be solved exactly in constant time by taking the smallest of the nine. 

Suppose A" has non-negative slope. Using the notation of the previous case and a 
similar argument, A > and 

< 9 , (26) 
arctan(A) < 9 , (27) 
9 < ir-u . (28) 



We obtain the same equations (24) and (25), and therefore the candidates are the same 
nine ones. So it can be solved exactly in constant time by taking the candidates that 
minimize n(9). 

Suppose A" is vertical. The general equation of A is x = s for s G 1. < u < |-7r 
otherwise there is no solution. Using the notation of the previous cases, there are three 
subcases to consider: (3.1) s < 0, (3.2) s G ]0, 2r sin(co>)[ or (3.3) s > 2rsin(o;). 

(3.1) Suppose A" is vertical and s < 0. Trivially < 9 < it — u, but there is one more 
restriction, v must be between q and b, hence 9 < — uj. It all sums up to 

< 9 , (29) 

9 < . (30) 

By elementary trigonometry, we get 

q = (2r cos(6 l ) sin(6 l + u), 2r sin(6 l ) sin(6 l + u)) 

cos(u;)cot(6') + sin (a; ) cot ^ 60 cos(co>) + sin(w)cot(6') 

ZT ZT 

cot 2 (#) + 1 ' cot\9) + 1 

Then, finding the general equation of the line through q and v , and of the line 
through q and w, we can calculate the general coordinates of b and c: 



cot(0) 

(s — 2r sin(ci;))(sin(a;)cot(^) + cos(a;)) 
cos (a;) cot (9) — sin(a;) 

We note X{9) = cot(#), therefore 

/ sin(a;)X 2 + cos(u)X sin(co>)X + cos(c<;' 
q = 2r — ,2r- 



6 = I. 



X 2 + 1 ' X 2 + 1 

s 



(s — 2r sin (a;) ) (sin (a;) X + cos(a;)) 
$ 

cos(c<;)X — sin(a;) 
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Here is the area a of Aqbc. 

a{X) 

1 



sinu 



((2rsin(w) - s)X 2 + 2rcos(w)X - a) 
X 2 (X 2 + 1) 



x 



((2r sin(w) - s)X 2 + 2r cos(u;)X - s) 2 
cos(w)X - sin(w)) 2 (X 2 + 1) 

((2r sin(w) - s)X 2 + 2r cos(w)X - s) 2 
1 |X(cos(u;)X - sm(u))\ (X 2 + 1) 
((2r sin(w) - s)X 2 + 2r cos(w)X - s) 2 
X(cos(a;)X - sin(w))(X 2 + 1) 



sinw 



sin (u 



The reason why the absolute value | • | disappeared is twofold. Firstly, by the 
hypothesis, < 9 < |7r, so X = cot(6 l ) > 0. Secondly, by (30), 9 < \ix — u, 
therefore 



9 < -7T - U) 



cot(6>) > cot y-n — cu 

cot (9) > tan(tu) 

cos(u;)cot(0) > sin (a;) 

cos(cu)cot(#) — sin (a;) > 

cos(a;)X — sin(w) > . 



Now we need to find for what values of X where cr(X) is minimum, wich means 
that we need to find for wich X does cr'(X) = 0. 

o'(X) 

= — sin(a;) ((2r sin(a;) — s)X 2 + 2r cos(a;)X — s) x 

((-s sin(w) + 2r sin 2 (w) + 4r cos 2 (w))X 4 - 2 cos(w)(s + 3r sin(w))X 3 

+6r sin 2 (o;)X 2 + 2 cos(o;)(r sin (a;) — s)X + ssin(a;)) -j- 
(X 2 (cos(w)X - sin(w)) 2 (X 2 + l) 2 ) 

Therefore cx'(X) = if and only if 

(2r sin(aj ) - s)X 2 + 2r cos(w)X - s = (31) 
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or 

(s sm(u) + 2r sm\oo) + 4r cos 2 (u))X 4 - 2 cos(w)(s + 3r sin(w))X 3 
+6r sm\u)X 2 + 2 cos(w) (r sin(w) - s)X + s em.(u) = . (32) 

The first equation is quadratic in X = cot(#), therefore it can be solved in constant 
time for X and then it remains to solve for 9. The second equation is quartic in 
X = cot(#), therefore it can be solved in constant time for X and then it remains 
to solve for 9. 

Actually, what we are really looking for is fi'(X(9)) = 0. Since [i'(X(9)) = 
n'(X)X'(9) = /i'(X)csc(6')cot(6 l ), 9 = \n is also a candidate. And of course, the 
extremities of the domain of 9 are also candidates. 

Overall, we get at most nine candidates for minimum (at most two from the 
quadratic equation, at most four from the quartic equation, 9 = |7T and the 
extremities.) So it can be solved exactly in constant time by taking the candidates 
that minimize /i(#). 

(3.2) Suppose A" is vertical and s G ]0, 2r sin(a;)[. In this case there is no solution. 

(3.3) Suppose A" is vertical and s > 2rsin(co>). Using the notation of the case 'A" is 
vertical and s < 0" and a similar argument, 

< , (33) 
9 < n -co . (34) 



We obtain the same equations (31) and (32), and therefore the candidates are 
the same nine ones. So it can be solved exactly in constant time by taking the 
candidates that minimize n(9). □ 



Lemma 12 describes how to compute the minimum enclosing triangle when q is forced to 



stay on one of the vertices of P and m is forced to stay on another one of the vertices of P. 

Lemma 12. Let q and Vj be points. It is possible to find the triangle Aqbc of minimum 
area such that Vj is the midpoint of be and Zbqc = uo in 0(1) time, provided that qb and qc 
are restrained to a closed wedge W strictly included in a 2u-wedge W such that the angle 



bisector ofW contains Vj (refer to Figure 22). 



PROOF: Without loss of generality, q = (0,0) and Vj = (//cot (a;), /x) with fi > 0. Let A 
(respectively A') be the line of general equation x— cot(9)y = (respectively x— cot(9+u)y = 



0) for < 9 < uo. (refer to Figure 23). The w-wedges created by A and A' for < 9 < to 
are all possible uo- wedges of vertex q that enclose Vj. Therefore, we have to find the value 
of 9 such that the triangle Aqbc (where b G A, c e A' and Vj is the midpoint of be) is of 
minimum area. 

Let A" be the line through q that is perpendicular to A. The general equation of A" is 
cot(9)x + y = 0. 
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Fix 9 and apply Lemma [5] where v := Vj and the coordinate system is the one created 
by A and A". Let s (respectively t) be the coordinate of Vj with respect to the A axis 
(respectively with respect to the A" axis). Then, following the proof of Lemma |5j let 
b (respectively c) be the point of A (respectively of A') such that \qb\ = 2(s — tcot(u;)) 
(respectively \qc\ = 2icsc(a>)). 

Here is the area a of Aqbc. 

= ^l9 6 ll9 c l sin M 

= 2^( s ~ tcot{u))2tcsc(uS) sin(cj) 

= 2t(s -tcot(w)) 

= 2(cos(6>) — sin(0)cot(u;))/i x 

((sin(6>) + cos (9) cot (w))}J> — (008(6*) — sin(6 l )cot(a;))/xcot(c<;)) 
= 2/i 2 csc 3 (w) sin(w - 9) sin(0) 

Now we need to find for what values of 9 is cr(^) minimum, wich means that we need to find 
for what 9 does o~'(8) = 0. 

a' (9) = 2/i 2 csc 3 (w) sin(w - 20) 

so <y'{9) = if and only if 



sin(w - 29) = 



therefor^] = | 



2 w. 



Is 9 = \u a minimum or a maximum? 



= -4/i 2 csc 3 (w) cos(w - 20) 



so 



a" I -u; 



1 
2 

2„„„3, 



= ^/Tcsc^u;) 
< 

because < u; < it. Hence 9 = \ti is the unique (local and global) maximum. 

Therefore, if the domain of 9 is restricted to < <pi < 9 < <p2 < tt, then the minimum of 
fi(9) is either at 9 = <pi or at 9 = <p 2 and it can be computed in constant time. □ 



Lemma 13 describes how to compute the minimum enclosing triangle when q is forced to 



stay on one of the vertices of P and m is forced to stay on one of the edges of P. 

Lemma 13. Let q be a point and A" be a line. The triangle Aqbc of minimum area such 



that b,c G A" and Z.bqc = u is isosceles and can be found in 0(1) time (refer to Figure 24) 



4 It corresponds to the case where Aqbc is isosceles. 



30 



Figure 24: Illustration of Lemma 
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PROOF: Without loss of generality, A" is horizontal, q is above A" and b is to the left of c. 
We need two intermediate results which follow from elementary geometry. 
Intermediate Results 

1. Let Aqbc and Aq'b'c' be two isosceles triangles such that \qb\ = \qc\, \q'b'\ = \q'c'\ and 



Zbqc = Zb'q'c' = u (refer to Figure 25). Denote by h q (respectively by h q /) the height 




Figure 25: Intermediate result 1. 

of Aqbc (respectively of Aq'b'c') relative to be (respectively to b 
\b'd\ > ' 



'd) 



If h. 



> h q , then 



2. Let Aqbc be an isosceles triangle such that \qb\ = \qc\ and Zbqc = u>. Let Aq'bc be 
a triangle such that Zbq'c = u (refer to Figure 26). Denote by h q (respectively by 



h q i) the height of Aqbc (respectively of Aq'bc) relative to be (respectively to be). Then 
h q > < h q . 

We are now ready for the main demonstration. 



6',c' G A", Zb'qc' = u and b' is to the left of d (refer to Figure 27) 
isosceles triangle such that \q'b'\ 



\q'c'\ 



and Zb'q'c' 



Let Aqb'd be any triangle such that 
Let Aq'b'c' be the 

LU. 



Denote by h q (respectively by 



h q i) the height of Aqbc (respectively of Aq'b'c') relative to be (respectively to b'd). By 
Intermediate result 2, h q > > h q , so \b'c'\ > \bc\ by Intermediate result 1. Since Aqb'd and 
Aqbc have the same height, then Aqbc is of minimum area. □ 
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Figure 26: Intermediate result 2. 




Step 4 For each pair of consecutive event points, compute the minimum triangle enclosing 
P and having a vertex on the corresponding arc of a circle. 

Since there are 0(n) event points, Step 4 takes 0(n). 

6 Putting it All Together 

The vertex that subtends the angle u of the optimal triangle enclosing P has to be on the 
w-cloud Vt. VL was computed in Step 1 in 0(n) time (refer to Section [2]). In Step 2 (refer 
to Section |3|, we fixed an w-wedge W and computed the minimum triangle that can be 
constructed with W in 0(n) time. This triangle is such that the midpoint m of its third side 
touches P (refer to Proposition [3] and Corollary [6]) . Then we divided Vt into a linear number 
of pieces in Step 3 (refer to Section Within each of these pieces, an optimal triangle can 
be computed in 0(1) time. This was done in Step 4 (refer to Section [5]). 

Step 5 Find the smallest triangles among those calculated in Step 4. 

Since there are 0(n) event points, Step 5 takes 0(n) time. Since each step takes no 
longer that 0(ri) time, then the algorithm takes 0(ri) time. If the input is a set of points, 
compute the convex hull and then apply this algorithm. In this situation, the computation 
clearly takes 0(n log(n)) time. 

We summarize the final algorithm. 
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Algorithm 14. (Minimum Enclosing Triangle with a Fixed Angle) 



• INPUT: A finite set S of points in the plane and an angle < 00 < n. 

• OUTPUT: All triangles with minimum area having an angle of u that enclose S. 

0. Compute the convex hull of S and denote it by P. We denote the edges and the vertices 
of P in clockwise order by and v,i for < i < n — 1 (all index manipulation is modulo 
n). 

1. Compute the co-cloud around P and denote it by (refer to Section^. Q consists of 
n' = 0(n) circular arcs that we denote in clockwise order by Tj for < i < n' — 1. The 
intersection point of Tj and r i+1 is denoted by for < i < n' — 1 

2. Let q G T be such that q = u$. Consider the u-wedge W = W(u, q, A, A') that encloses 
P. Apply Algorithm^ with P and W (refer to Section^). Let Aqbc be the ouput of 
Algorithm^ and denote by m the midpoint of segment be. 

3. Move q clockwise along Q and maintain W, b, c and m as defined in^ Collect all of 
the following three types of event points (see Section^ for formal definition) 

Type 1 q is on the intersection point of two consecutive circular arcs of Q for an i 
with < i < n' — 1 . 

Type 2 q is such that the third side be of the triangle is on an edge e« of P (for an i 
with < i < n — 1) and the midpoint m of be is on the first vertex Vi of ( when 
the vertices of P are considered in clockwise order). 

Type 3 q is such that the third side be of the triangle is on an edge e, of P (for an 
i with < i < n — 1) and the midpoint m of be is on the last vertex Vi+\ of e, 
(when the vertices of P are considered in clockwise order). 

4- For each pair of consecutive event points computed in^ which we denote by (q',q"), 
there is a single arc of a circle. This circular arc might be reduced to a single point if 
two consecutive event points are on one of the vertices of P. When q moves on such 
an arc, either m is forced to stay on one of the vertices of P or m is forced to stay on 
one of the edges of P (see Sections^ and^for complete discussion). For each pair 
(q', q") of consecutive event points, 

• if the circular arc between q' and q" is not reduced to a single point and if m is 



forced to stay on one of the vertices of P, store the triangle defined by Lemma 10 
(refer to Section^. 

If the circular arc between q' and q" is not reduced to a single point and if m is 



forced to stay on one of the edges of P, store the triangle defined by Lemma \Tl 
(refer to Section^. 
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If the circular arc between q' and q" is reduced to a single point and if m is forced 



to stay on one of the vertices of P, store the triangle defined by Lemma 12 (refer 
to Section^. 

If the circular arc between q' and q" is reduced to a single point and if m is forced 



to stay on one of the edges of P, store the triangle defined by Lemma 13 (refer to 
Section^. 

5. Return the smallest triangles among those stored in^\ 



7 A Note on the Complexity of the Solution 

Notice that one of the cases in Section [5] required us to find the roots of a fourth degree 
polynomial (refer to Lemma 11 ). One may ask whether or not this is necessary or if one can 



avoid finding the roots of such a polynomial to solve this problem. In this section, we show 
that it is unavoidable in certain situations, by providing an example of a set of points where 
the optimal solution lies on the root of an irreducible fourth degree polynomial. 

Consider the quadrilateral [abed] where a = (0,0), b = (2,0), c = (2,— |) and d = 

^— 4 4v ^T x , 4 v/ ^~ 12 j , and take u> = \n (hence we are looking for the minimum enclosing right 

triangle). It turns out that the optimal right triangle is such that the right angle is on T 



(refer to Figure 28) and the hypotenuse is on cd. Therefore, we need to solve (25) with 



Optimal triangle 




Figure 28: Example where the optimal solution involves the roots of a quartic equation. 

[i := —1, A := — \ and r := 1. It leads to the following quartic equation. 

p A (X) = 13X 4 - 92X 3 + 45X 2 + 12X - 62 = 
This equation admits two real solutions Xi and X 2 . Note 



a = V -1107 + 6^34134 
34 



1105 221 

1726 + a 

a 3 



then 



X 2 



23 
13 
23 
13 



26 
1 



P 



26^+2 



3452 
169 


17 

a - 

39 


85 
13a 


+ 


136796 
169/3 


3452 


17 

a - 

39 


85 


+ 


136796 


169 


13a 


169/3 



-0.761694 



6.543373 



Solving cot (9) = X\ and cot(0) = X%, then testing against a{9) (refer to the proof of 
Lemma 11), one finds that 9 ~ arccot(6. 543373) ~ 0.048273 ir minimizes u{9). 



Since Pa(X) and its resolvent cubic 

p{X) = X 3 - 



45 
13* 



X 2 



2120 
l69~ 



X 



377816 
2197 



are irreducible over Q [X] , the algebraic expressions for the roots of p4 must be written with 
square roots and cubic roots. Moreover, they cannot be simplified (see |DF03j ). 

Therefore, in general, we cannot avoid square roots nor cubic roots in the computation 
of the minimum enclosing area triangle with a fixed angle. 



8 Conclusion 

We have shown how to compute all triangles of minimum area with a fixed angle < u < it 
that enclose a point set. It would be interesting to see if some of these techniques generalize 
to other settings or other optimization criteria. For example, finding the smallest area 
tetrahedron with a fixed solid angle of a set of points in three dimensions. 
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